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We study the vortex distribution of the wave functions minimizing the Gross Pitaevskii energy for 
a fast rotating condensate in the Lowest Landau Level (LLL): we prove that the minimizer cannot 
have a finite number of zeroes thus the lattice is infinite, but not uniform. This uses the explicit 
expression of the projector onto the LLL. We also show that any slow varying envelope function 
can be approximated in the LLL by distorting the lattice. This is used in particular to approximate 
the inverted parabola and understand the role of "invisible" vortices: the distortion of the lattice is 
very small in the Thomas Fermi region but quite large outside, where the "invisible" vortices lie. 

PACS numbers: 03.75.Hh, 05.30.Jp, 67.40.Db, 74.25.Qt 



The fast rotating regime for a Bose Einstein condensate 
in a harmonic trap, observed experimentally in 
displays analogies with type II superconductors behav- 
iors and Quantum Hall Physics. However, some different 
features have emerged and are of interest, in particular 
due to the existence of a potential trapping the atoms. 

A quantum fluid described by a macroscopic wave func- 
tion rotates through the nucleation of quantized vortices 
0. For a condensate confined in a harmonic potential 
with cylindrical symmetry around the rotation axis, a 
limiting regime occurs when the rotational frequency fl 
approaches the transverse trapping frequency: the cen- 
trifugal force nearly balances the trapping force so that 
the size of the condensate increases and the number of 
vortices diverges. The visible vortices arrange themselves 
in a triangular Abrikosov lattice. The system is strongly 
confined along the axis of rotation, and it is customary to 
restrict to a two dimensional analysis in the x — y plane. 
We will call z — x + iy. The hamiltonian is similar to 
that for a charged particle in a magnetic field: for rota- 
tional angular velocities just below the transverse trap 
frequency, the wave function of the condensate can be 
described using only components in the Lowest Landau 
Level (LLL): 
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where $o is a normalization factor and the z, are the 
location of the vortices. In rescaled units, the reduced 
energy in the LLL is H 0, Q 

£LLLW = J[(l~n)\z\^ + ^}d 2 r (2) 

under Jd 2 r\^\ 2 = 1, where f2 is the rotational velocity, 
the transverse trap frequency is scaled to 1, and G mod- 
els the interaction term: G — Ng/(dV2w), where g is the 



two body interaction strength and d is the characteris- 
tic size of the harmonic oscillator in the direction of the 
rotation. 

In the absence of a confining potential, the problem is 
reduced to the one studied by Abrikosov |8( for a type 
II superconductor and the minimizer is a wave function 
with a uniform triangular lattice 0; its modulus van- 
ishes once in each cell and is periodic over the lattice. 
The presence of the confining potential is at the origin of 
a slow varying density profile, which can be described as 
the mean of the modulus of the wave function on many 
cells. Ho p] predicted that for a uniform lattice, the 
smoothed density profile is a gaussian. Various contri- 
butions [a, IjJ llfl ] then pointed out that the energy can 
be lowered if this smoothed density distribution is an in- 
verted parabola rather then a gaussian. This type of den- 
sity profile can be achieved either by taking wave func- 
tions with a uniform lattice but with components outside 
the LLL [a] or by remaining in the LLL and distorting 
the lattice. The study of the distortion has been the fo- 
cus of recent papers @, 0, 0] and raises the issue of 
the optimal vortex distribution. In the LLL description, 
there are two kinds of vortices : the "visible vortices", 
which lie in the region where the wave function is sig- 
nificant (for instance inside the Thomas Fermi region in 
the case of the inverted parabola) , and the "invisible vor- 
tices" which are in the region where the modulus of the 
wave function is small. The visible vortices form a reg- 
ular triangular lattice, while the invisible ones seem to 
have a strong distorted shape, whose distribution is es- 
sential to recreate the inverted parabola profile inside the 
LLL approximation. These latter are not within reach of 
ex per imental evidence, but can be computed numerically 
[rllllj. An important theoretical question is the distribu- 
tion of these invisible vortices, their number, or an esti- 
mate of how many of them are necessary to approximate 
the inverted parabola properly inside the LLL. 
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Our main result is to prove that in order to minimize 
the energy in the LLL, there is a need for an infinite 
number of vortices. The main tool that we use is an 
explicit expression of the projector onto the LLL. This 
projector also allows us to approximate any slow varying 
density profile by LLL wave functions. 

Projection onto the LLL and infinite number of 
zeroes We define a small parameter e = yl — and 
make the change of variables ip(z) — y/e^ly/ez), so that 
the condensate is of size of order 1 and the lattice spacing 
is expected to be of order yfe. The energy gets rescaled 
as £ L ll(^) = eE LLL (ip) where 
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d 2 r. 
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Moreover, if) belongs to the LLL so that f(z) — 
ip(z)e\ z \ / 2e is a holomorphic function and thus belongs 
to the so called Fock-Bargman space 



T = j/ is holomorphic , J |/| 2 e |z|2/e d 2 r < ocj 



(4) 



Let us point out that such a function / is not only de- 
termined by its zeroes and normalization factor as in Q , 
but also by a globally defined phase, which is a holomor- 
phic function. The space J- is a, Hilbert space endowed 
with the scalar product (f,g) = J f(z)g(z)e~^ z ^'^ E d 2 r. 
The point of considering this space is that the projection 
of a general function q(z, z) onto T is explicit, and called 
the Szego projector [12|,LUj : 



n(s) = - 

ne 



- g{z' , z')d 2 r' 



(5) 



If g is a holomorphic function, then an integration by 
part yields II (g) — g. 

If one considers the minimization of Elll(i[)) without 
the holomorphic constraint on /, then the minimization 
process yields that |z| 2 + G|?/>| 2 — /i = 0, where fi is the 
chemical potential due to the constraint J \if>\ 2 — I, so 
that \ip\ is the inverted parabola 
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(6) 

The restriction to the LLL prevents from achieving this 
specific inverted parabola since ^e' 2 ' l 2e cannot be a 
holomorphic function. The advantage of the explicit for- 
mulation of the projector II is that it allows us to derive 
an equation satisfied by if) or rather / when minimizing 
the energy in the LLL. A proper distribution of zeroes 
can approximate an inverted parabola profile but is go- 
ing to modify the radius R by a coefficient 6 1 / 4 coming 
from the contribution of the vortex lattice to the energy. 



Theorem 1 If f G T minimizes 
E(f) = 



d 2 r 
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under J \f\ 2 e ' z ' l £ d 2 r — 1, then f is a solution of the 
following equation 



n((M 2 + G|/| 2 



M 2 /e 



-t*)f =o 



(8) 



where /i is the chemical potential coming from the mass 
constraint. 

Note that given the relation between / and if>, E(f) and 
Elll(iP) are identical. Equation JHJ comes from the fact 
that for any g in T with (/, g) — 0, if / minimizes E, 
then we have 



2 ff/e^ |2/e + f|/| 2 ff/ e - 2|z|2/£ 



d 2 r = 



(9) 



and we use the scalar product in T and the definition of 
the projector to conclude. 

The equation for the minimizer allows us to derive that 
this minimizer cannot be a polynomial: 

Theorem 2 If f G T minimizes E, then f has an infi- 
nite number of zeroes. 

We are going to argue by contradiction and assume that 
/ is a polynomial. 

1. The proof first requires another formulation of JSJ. 
The projector II has many properties m partic- 

ular, one can check, using an integration by part in the 
expression of II, that II £ (|z| 2 /) = zed z f + ef. As for the 
middle term in the equation, one can compute that if / is 

a polynomial, U £ (e"^ \f\ 2 /) = n e (V^ |/| 2 ) II E / 

= n £ (7(z))rL(e-^/ 2 ) = f{ed z )U e (e--f f 2 ). A 
simple change of variable yields IJ e ( e ^/ 2 ) ( z ) = 

(ire)- 2 fe- z ~~ 2 ^' 12 f 2 (z')d 2 r' 

= ^ (^(72)) ^ = y 2 ^- xhus ' we find the foi " 

lowing simplification of (jHJ: 
G 



zed z f + -f(ed z )[f 2 (z/2)] - (jjt - e)f = 0. 



(10) 



2. Now we assume that / is polynomial of degree n 
and a solution of l |10| l. We are going to show that there 
is a contradiction due to the term of highest degree in 
the equation. Indeed, if / is a polynomial of degree n, 
then (ed z ) k [f 2 (z/2)} is of degree In - k. But (JTHJ im- 
plies that f(ed z )[f 2 (z/2)] is of degree n, hence / must 
be equal to cz n . This is indeed a solution of (|f Oj l if 
ne + G\c\ 2 £ n {2n)\/{2 2n+1 n\) - fi + e = 0. Using that 
1, we find that \c\ 2 Tie n+1 n\ = 1. The 
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Stirling formula provides the existence of a constant cq 



3 



such that ne + coG/(2Tte^/n) < /i. For the minimizer, fi is 
of the same order as the energy, thus of order 1, so that if 
e is too small, no n can satisfy this last identity hence the 
minimizer is not a polynomial. A similar argument can 
be used to check that, if / is more generally a holomor- 
phic function in J 7 , then it cannot have a finite number 
of zeroes. The detailed proof will be given in 

El- 
Approximation of a slow varying profile by the 

LLL The Abrikosov solution The Abrikosov problem [j| 
consists in minimizing the ratio (|w| 4 ) / (|w| 2 ) over peri- 
odic functions, where (.) denotes the average value over a 
cell, for functions u obtained as limits of LLL functions. 
The minimum is achieved for u = u £ 0, e 2 ^/ 3 ) where 



U £ (z, T) = e -|*l 2 /2e /£ ( Zj T) /e( t) = e *'/2*e( / ^ T ) 

V 7T£ 



(11) 

and for any complex number r = tji + itj , 

9( W;T ) = _ J2 (-l) n e i7rr( " +1/2)2 e (2n+1)7rra . (12) 

n— — oo 

The O function has the following properties 

e(v + k + It, t) = {-l) k+l e~ 2Mv e- MT Q{v, r) (13) 



so that |u e (z,t)| is periodic over the lattice y 77Z 

^J~^Zt, and vanishes at each point of the lattice. With- 
out loss of generality, one can restrict r to vary in |t| > 1, 
— 1/2 < tr < 1/2: this is equivalent to require that the 
smallest period for is 1 and along the x axis (see H3) 
and any lattice in the plane, can be obtained from one of 
these by similarity. 

For any r, f s given by Ijlll) is a solution of 

n(|/e| 2 e " |z|2/£ /s) = Xrfe, With \ T = (M 2 ) b(T), (14) 
and 

(kl 4 ) 



b(r) 



^2 e _7r|fer 



(15) 
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This expression can be obtained using arguments in 
[lfij |. The minimal value of b(r) ~ 1.16 is achieved for 
t = e 2l7r / 3 , that is for the triangular lattice Q: in 0], it 
is argued that one can restrict to tr = —1/2, and vary 77 
in (1/2, y/3/2). Accepting this restriction, they compute 
the variations of b which depends on a single parameter 
and is indeed minimal for the triangular lattice. In [l7j . 
we prove that this restriction is rigorous using the de- 
scription of these lattices by varying r for |t| = 1 and 

TR e (-1/2,0). 

If one compares (|14f) and (JSJ, one notices that they 
only differ by the term II(|z| 2 /) = ezd z f + ef, which 



is negligible on the lattice size, but plays a role on the 
shape of the density profile. 

The role of the confining potential A natural candidate 
to approximate any slow varying profile a(z, z) is to take 
a(z, z)u e (z, t), where u £ is the periodic function defined 
in l|lll) . Of course, such a function is not in the LLL, 
but can be well approximated in the LLL by / Q e~' z ' / 2e 
where f a = II(a/ e ), II is the projector onto the LLL (0) 
and f £ comes from (|11H . Estimating the energy of /" 



yields £(r)-/[M 2 H 2 (k| 2 ) 



Gb(r 



Ce 1 / 4 . This computation uses calculus on II |l7|. and 
that u e and a do not vary on the same scale, hence the 
integrals can be decoupled. The contribution of u £ to the 
energy is through the coefficient 6(r), which is minimum 
for t = e 2i7r / 3 . 

Using pseudo differential calculus, one can show [l7| . 
when e is small, that f a is very close to au e : the error is 
at most like e 1 ' 4 if a is not more singular than an inverted 
parabola. In particular, when a is an inverted parabola, 
this implies that in the Thomas Fermi region, the distri- 
bution of visible vortices is almost that of the triangular 
lattice since au s is a good approximation. Outside the 
support of the inverted parabola, where f a is very small, 
one can check that the density of distribution of zeroes 
of f a decreases like 1 / \z I for large \z\. Contrary to what 
was explained in 0, [la | , it is not a small distortion of 
the lattice which results in large changes in the density 
distribution, but a very specific and far from uniform 
distribution of the invisible vortices (outside the Thomas 
Fermi region) which allows to approximate an inverted 
parabola. 

The special shape of the inverted parabola comes out if 
one wants to approximate the equation of the minimizer 
of the energy: for any A, we can prove that 

n((|z| 2 + G|n 2 e - |2|2/e - A)/ Q ) + Ce 1/4 
~ n((|z| 2 + Gb(r) (M 2 ) H 2 - X)af e ) (16) 



where C only depends on bounds on a. In other words, 
in the equation for / Q , one can separate in the term 
|ya|2 e -|z| /e contributions due to the lattice and to 
the profile. The right hand side of (|16|) is zero if a is the 
inverted parabola 



irRl (\u £ 



1 - 



R 2 



i?,n — 



2Gb(r) 



and A = i?g, so that f a is almost a solution of ©, up to 
an error in e 1 / 4 . 

Variations of the lattice This approach can be used 
to study the variations in energy due to deformations 
of the lattice. The triangular lattice, corresponding to 
t 1 = e 2l7r / 3 , is such that the Hessian of 6(t) is isotropic 
(~ 0.68/d). Two lattices close to each other can be de- 
scribed by two close complex numbers r 1 and r 2 ; the dif- 
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ference in energy between E(f a (.,T 1 )) and E(f a (.,T 2 )) 
is at leading order 



G d 2 b 



2 1 2 



|a| 4 (| U£ | 2 ) 2 d 2 r 



0.68G 



2|2 



This computation justifies the approach which consists 
in decoupling the lattice contribution from the profile 
contribution in the energy [lij but, given the definition of 
f a using IT, it relies on strong deformations of the lattice 
for points far away from the Thomas Fermi region. For 
a shear deformation for which Uij are the components of 
the deformation tensor, r 2 — r 1 = i\[iu xy . The elastic 
coefficient C2 is defined by the fact that the difference 
in energy should be 4C2U 2 y . This separation of scales 
allows to compute C 2 ~ 0.68G/(4tt j R 2 ) (see also 0E3) 
and relate it in BEC to the same one computed for the 
Abrikosov solution. 

Approximation by polynomials and modes An 
interesting issue, especially for the computations of 
modes, is to get an estimate of the degree of the poly- 
nomial which could approximate / Q , since this function 
has an infinite number of zeroes. We can prove that 
as the degree of the polynomial gets large, the minimum 
of the energy for the problem restricted to polynomials 
(and the computation of modes) is a good approxima- 
tion of the full problem. The convergence rates that we 
obtain are not satisfactory yet. We believe that a good 
degree should be k/e, with k > Rq and Rq is the radius 
of the inverted parabola. Given that the volume of the 
cell is 7T£, k = Rq would correspond to having only the 
visible vortices. Numerical simulations indicate that a 
sufficient number of invisible vortices is needed to recre- 
ate the inverted parabola profile Q • There are two types 
of invisible vortices: those close to the boundary of the 
inverted parabola which contribute to the bulk modes 
and those sufficiently far away which produce single par- 
ticle excitations as explained in An open issue is to 
understand the location of these latter invisible vortices; 
some simulations suggest that they lie on concentric cir- 
cles, but then the density of these circles should be very 
low to match our predicted global vortex density far away 
which behaves like 1/|^|- We have performed numerical 
simulations with 17 = 0.999 and G — 3: this fixes the 
number of visible vortices to 30, and we vary the number 
of total vortices N . One needs at least TV = 52 (that is 22 
invisible vortices) to properly approximate the inverted 
parabola, the energy minimizer and the bulk modes. The 
distortion of the lattice is small at the edges but large at 
large distances. For N too small, some modes do not ap- 
pear (see Figure while for N very large, one expects 
higher modes that 0, [2(| interpret as single particle 
modes. 

Conclusion We have shown that for the minimizer 
of the Gross Pitaevskii energy in the LLL, the lattice of 
vortices is infinite, but not uniform. Any slow varying 
profile can be approximated in the LLL by distorting the 




FIG. 1: Number of modes n/ntot having a lower energy than 
a given energy e for G — 3, Q = 0.999 



lattice. This is proved using an explicit expression for the 
projection onto the LLL. Our results also give an insight 
on the elastic coefficient C2 and the approximation of the 
minimizer and modes by polynomials. 
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